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Galaxy mergers at high resolution: 

From elliptical galaxies to tidal dwarfs and globular clusters 
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Numerical simulations of galaxy mergers are a powerful tool to study these fundamental events in the hierarchical built- 
up of galaxies. Recent progress have been made owing to improved modeling, increased resolution and large statistical 
samples. We present here the highest-resolution models of mergers performed so far. The formation of a variety of sub- 
structures ranging from kinematically decoupled cores to globular-like clusters is directly resolved. In a resolution study, 
we show that the large-scale structure of elliptical-like merger remnants can be affected by the resolution, and a too mod- 
est resolution may affect the numerical predictions on the properties of major merger remnants: understanding precisely 
which kind of event or succession of events has formed the various types of elliptical galaxies remains an open challenge. 
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1 Introduction 

For more than two decades, numerical simulations have 
been a powerful tool to study the dynamics of interact- 
ing galaxies and the properties of galaxy merger remnants. 
Some of the main results obtained from numerical models 
range from the formation of elliptical-like galaxies with re- 
alistic r^/^ density profiles in major mergers of massive spi- 
rals (Barnes 1992, Naab & Burkert 2003, Bournaud et al. 
2005), the effect of minor mergers of spiral galaxies with 
small companions (e.g., Walker et al. 1996), and the trig- 
gering of star-formation by galaxy interactions and mergers 
(Mihos & Hernquist 1994, Cox et al. 2006, di Matteo et al. 
2007). 

Recent progress have been made by studying galaxy in- 
teractions in a wider range of physical conditions includ- 
ing more complex (and more realistic) cases like: the (re- 
)merging of early-type galaxies that are themselves rem- 
nants of past major mergers (Naab et al. 2006), sequences of 
repeated minor mergers (Bournaud et al. 2007a), or merg- 
ing galaxies embedded in a larger- scale group or cluster 
environment, showing a significant influence of the ram- 
pressure and large-scale gravity field on interacting galaxies 
(Kapferer et al. 2008, Martig & Bournaud 2008). 

We present here recent progress in terms of numeri- 
cal resolution, in simulations of "standard" mergers of two 
massive spiral galaxies merging together, with no multiple 
mergers and no large-scale environment accounted for (see 
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Section 2). While the remnants of major mergers in nu- 
merical simulations do globally resemble elliptical galax- 
ies, whether the major merger scenario can really explain 
the exact properties of most real elliptical galaxies or not is 
still an open question (see Burkert et al. 2007). This raises 
the question of whether or not most simulations of galaxy 
mergers have a high-enough resolution to predict accurately 
the detailed properties of merger remnants: this is studied in 
Section 3, where we show that a too modest resolution may 
significantly affect the properties of numerical merger rem- 
nants. Our recent very-high resolution simulations also en- 
able to study the smallest structures formed in galaxy merg- 
ers, like decoupled cores, tidal dwarfs, and super star clus- 
ters: results on these aspects will be presented in Section 4. 



2 The simulations and resolution 

The simulations presented here were performed with a 
particle-mesh code based on a FFT Poisson solver for the 
gravity. The gas dynamics is modeled with a sticky-particle 
scheme, and a Schmidt law is used for star formation. A 
thorough description of the code and initial conditions can 
be found in Bournaud et al. (2008) and references therein. 
The stellar disks have initial masses of 2 x 10^^ M© and 
we use Burkert (core) profiles for the dark halos. In the fol- 
lowing we present simulations obtained at three different 
resolutions, all with the same 0.5 Gyr timestep: 

- the low resolution has a softening of 180 pc (up to 
25 kpc from each galaxy center, and degraded at larger 
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Fig. 1 Projected stellar density maps for a major merger 
between two spiral galaxies. The low, medium and high res- 
olution models of the same merger are viewed with the 
same projection angle. These snapshots show an area of 
16x16 kpc. 



distances) and uses 10^ particles/galaxy/component (the 
components are: stars, gas, and dark matter; this implies 
a total of 6 X 10^ particles). 

- the medium resolution has a softening of 80 pc and 10^ 
particles/galaxy/component (total of 6 x 10^ particles). 

- the high resolution has a softening of 32 pc and uses 
6 X 10^ particles/galaxy/component (total of 36 x 10^ 
particles). Too our knowledge, this is the highest resolu- 
tion ever achieved so far in this kind of simulations (see 
comparison with other works in Bournaud et al. 2008). 

3 Global properties of elliptical galaxies: a 
resolution study 

To test the standard proposal that ellipticals have been 
formed by galaxy mergers, and explore in more detail which 
sort of merger has formed them (binary major mergers, re- 
peated minor or major merger, wet/dry mergers, etc.), re- 
cent works have measured in detail the properties of such 
numerical merger remnants (references above). This raises 
the question of whether or not numerical simulations with 
typically a few 10^ particles per galaxy have converged in 
terms of the detailed properties of these "ellipticals". 

Here we present a wet 1:1 merger model: the two spiral 
progenitors contain 83% of stars and 17% of gas (plus a 
dark halo with a Burkert profile), and they have the same 
mass. The simulation was run at the low, medium and high 
resolutions defined above. The morphology of the relaxed 
merger remnant under the very same projection is shown for 
each resolution on Figure 1 . Among the many quantitative 
parameters measured (which will all be presented in Bois 
et al. 2009, in preparation), we show on Figure 2 the radial 
variations of the robust tracer of angular momentum A^^, 
defined by Emsellem et al. (2007), measured for a set of 
isotropically distributed projections. 

Our main conclusions follow, and we made sure (Bois et 
al. 2009) that these effects are not artifacts caused by some- 
what different initial conditions (defining the very same ini- 
tial conditions at different resolution is not trivial), do not 
affect only the gas dynamics and star formation but also the 
stellar and dark matter dynamics (sizeable results are found 
for dry mergers), and do not result simply from a change 
in the orbit (because the dynamical friction would not be 



equally resolved in the three cases), but really result from 
the way the violent relaxation processes are resolved: 

- At the three resolutions, the merger remnant resemble 
an elliptical galaxy, with about the same flattening (axis 
ratio) except in the central 3 kpc, and a similar r^/^-like 
density profile. The effective (half-mass) radius is 2- 
2.5 kpc, depending on the projection angle, in the three 
cases. 

- In more detail, significant variations can be seen in the 
isophotal shape (Fig. 1). The low resolution case is over- 
all disky with a major (~ 70°) isophotal twist. The 
medium and high resolution cases are boxy under the 
same projection, with little isophotal twist. 

- Important kinematical variations are seen (Fig. 2), the 
angular momentum profile of the low resolution case 
differs dramatically from the medium and high resolu- 
tion ones - even the two later cases differ from each 
other, but in a more modest way and the shape of their 
\r profiles is more comparable. In particular, within 
their effective radius, the low resolution model would 
be classified as a fast-rotator (as defined by Emsellem et 
al. 2007) while the medium and high ones appear to be 
slow-rotators: this is a major difference as this typically 
defines two different classes of early-type galaxies (Em- 
sellem et al. 2007), and the origin of the slowly-rotating 
ellipticals is a particular challenge for merger scenarios 
(as pointed out by Naab & Ostriker 2007). 

Overall, our low resolution model (180 pc, 10^ parti- 
cles/galaxy /component), which is not so "low" compared to 
recent works in the field, gives only a gross prediction of the 
properties of the merger remnant, and the merger/relaxation 
processes are insufficiently resolved to accurately predict 
the detailed properties of this "elliptical". Our medium reso- 
lution (80 pc, 10^ part/gal/comp), which is quite high com- 
pared to several recent works, gives more accurate predic- 
tions of the morphology and kinematics. It generally com- 
pares favorably to the high resolution case (32 pc, 6 x 10^ 
part/gal/comp - the highest to date), but some results may 
still vary with resolution, in particular the advanced param- 
eters used to deeply explore the structure of ellipticals like 
\r. This shows that resolution issues are a significant con- 
cern in attempts to compare in detail the properties of ob- 
served ellipticals with prediction of merger models, even in 
recent simulations with a rather high number of both stellar 
and gaseous particles. 

4 Structure formation in galaxy mergers: 
from decoupled cores to tidal dwarfs and 
globular clusters 

While the total stellar mass density of merger remnants is 
dominated by smooth spheroids (Fig. 1), interstellar gas fu- 
els the formation of many substructures of gas and young 
stars. Kinematically decoupled cores (KDCs) are frequently 
seen in our medium- and high-resolution simulations of 
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Fig. 2 Radial profiles of the angular momentum tracer Xr of a relaxed merger remnant, for the low, medium and high 
resolution models respectively. The many curves in each plot correspond to different projection angle, that were chosen to 
be isotropically distributed. 





Fig. 3 Stellar velocity field in a high-resolution wet 1 : 1 
merger remnant, showing a counter-rotating KCD. The 
black contours show the projected stellar mass density (log 
spacing), this snapshot shows a 10x12 kpc area. 

galaxy mergers (see Fig. 3). They form when a merger- 
driven gas inflow triggers the formation of a dense con- 
centration of young stars inside, typically, the central kpc. 
Such stellar systems formed during/after the merger largely 
contribute to the total stellar kinematics in the central re- 
gions, while the larger-scale system is dominated by stellar 
populations older than the merger. Torques that act on the 
gas structures during the galaxy collision frequently cause 
the central component to have a kinematic axis misaligned 
with the older stellar system, and the KDC can even some- 
times be counter-rotating w.r.t. the rest of the elliptical. Al- 
though our statistics are limited for the moment, we do not 
find KDCs to be preferentially associated to morphological 
decoupling like large isophotal twists. Note that major wet 
mergers are not necessarily the only mechanism producing 
kinematically decoupled components in early-type galaxies 
(Naab et al. 2007, di Matteo et al. 2008). 

The outskirts of merger remnants also contain a variety 
of young stellar substructures. Figure 4 displays the density 
map of the "young" stars formed during/after the merger, 
in a relaxed high-resolution merger remnant (see also Bour- 
naud et al. 2008) : one can notice streams and shells, and 
more than one hundred of small stellar objects. The three 
most massive stellar objects have masses of 10^~^ radii 
of about a kpc and are supported by rotation {V > cr - 
see objects A and B on Fig. 5): they resemble Tidal Dwarf 
Galaxies (TDGs). The many other, smaller objects have 
lower masses, ranging from 10^ (the resolution limit) 



Fig. 4 Projected mass density of the "young" stars, i.e. 
those formed during or after the merger, in a wet 1 : 1 high- 
resolution merger remnant. 

to a few 10'' M0. They are very compact with radii from 
~30 pc (the softening limit) to ~150 pc, and are tightly 
bound {Eg/2Ek < — 1, see Bournaud et al. 2008), sup- 
ported by random motions with little rotation (Y/cj < 1, 
see Fig. 5). These Super Star Clusters (SSCs) are thus likely 
progenitors of Globular Clusters (GCs). The formation of 
such SSCs/GCs in major mergers had been predicted for in- 
stance by Li et al. (2004) and is here directly resolved. 

TDGs are SSCs/GCs appear to be two different kind of 
objects: TDGs are large, rotating, found in the outer tidal 
tails. Their formation occurs when a giant molecular cloud 
or a dense region of the progenitor spiral disk is moved into 
a tidal tail by the interaction, following the mechanism pro- 
posed by Elmegreen et al. (1993). With large dark haloes, 
large regions of the outer disk can even be displaced into the 
outer regions of a tidal tail without being disrupted, and can 
form a new self-gravitating objects here (Due et al. 2004). 

The SSCs/GCs are formed in the outer tidal tails and 
other dense streams of gas at smaller radii, they have a 
much lower specific angular momentum and are more com- 
pact. This, together with a dip in the mass function between 
SSCs/GCs and TDGs (Fig. 6), implies that these are two dif- 
ferent categories of objects, and the three large TDGs in our 
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Fig. 5 Line-of- sight velocity and dispersion fields for sev- 
eral young stellar objects formed in the mergers: two TDGs 
(A, B) and three SSCs/GCs (C, D, E). 



log(M) 

Fig. 6 Mass function of the young bound stellar objects 
(SSCs/GCs and TDGs) formed in the major merger dis- 
played on Fig. 4. A significant dip separates SSCs/GCs and 
the three massive TDGs. The dashed histogram corresponds 
to a test simulation with modified gas dissipation parameters 
(see Bournaud et al. 2008). 

simulation are not just the high-mass end of the SSC distri- 
bution. The SSCs are formed by the fragmentation of dense 
gaseous structures, following the mechanism studied in par- 
ticular by Wetzstein et al. (2007). In that vein, the "tidal 
dwarfs" in Wetzstein et al. models do not resemble the three 
massive TDGs in our model, but rather to the high-mass end 
of the more compact SSCs. This does not mean that they 
should never be considered as tidal dwarfs: indeed, with 
masses of a few 10'' M©, they are really as massive as some 
dwarf galaxies. There could thus be two types of tidal dwarf 



galaxies: the high-mass end of SSCs/GCs, which would be 
compact spheroidal galaxies, and the larger, more massive 
and rotating TDGs, formed by a different mechanism and 
resembling the kpc-sized rotating TDGs observed in several 
interacting systems like NGC 5291 (Bournaud et al. 2007b). 
TDGs of this later type may still evolve into smaller dwarf 
spheroidals over a few Gyrs (Metz & Kroupa 2007). 

5 Conclusions 

High spatial and mass resolution enables the formation of 
small objects to be resolved in numerical models of galaxy 
mergers. In recent simulations, even the internal properties 
of small objects can be resolved - like the rotation curves of 
tidal dwarfs (see Bournaud et al. 2007b). The results of such 
high-resolution models is overall in agreement with theoret- 
ical expectations, in particular regarding the formation of a 
population of GCs during major wet mergers. 

High resolution is also a must to study accurately the 
large-scale structure of spheroidal galaxies formed in major 
mergers. Our recent resolution study shows that even sim- 
ulations with a total of one million particles and a spatial 
resolution of a hundred of parsec can still be significantly 
affected by the resolution. Whether major mergers can ex- 
plain the detailed properties of elliptical galaxies (boxyness, 
angular momentum distribution, anisotropy, etc..) is still 
largely unsolved, and understanding which type of merger 
would have formed the various kind of ellipticals remains 
an important challenge. 
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